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Abstract Recent analyses of nuclear decay data show evidence of variations 
suggestive of a solar influence. Analyses of datasets acquired at the Brookhaven 
National Laboratory (BNL) and at the Physikalisch-Technischc Bundcsanstalt 
(PTE) both show evidence of an annual periodicity and of periodicities with 
sidereal frequencies in the neighborhood of 12.25 year~^ (at a significance level 
that we have estimated to be 10""'^^). It is notable that this implied rotation 
rate is lower than that attributed to the solar radiative zone, suggestive of a 
slowly rotating solar core. This leads us to hypothesize that there may be an 
"inner tachocline" separating the core from the radiative zone, analogous to the 
"outer tachocline" that separates the radiative zone from the convection zone. 
The Ricger periodicity (which has a period of about 154 days, corresponding 
to a frequency of 2.37 year~^) may be attributed to an r-mode oscillation with 
spherical-harmonic indices I = 3,m = 1, located in the outer tachocline. This 
suggests that we may test the hypothesis of a solar influence on nuclear decay 
rates by searching BNL and PTB data for evidence of a "Rieger-like" r-mode 
oscillation, with Z = 3,m = 1, in the inner tachocline. The appropriate search 
band for such an oscillation is estimated to be 2.00-2.28 year~^. We find, in both 
datasets, strong evidence of a periodicity at 2.11 year~^. We estimate that the 
probability of obtaining these results by chance is 10^^^. 
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1. Introduction 

IJenkins et all (^005)) . iFischbach et ZlipOPg)) . and lJavorsek II et a?.l(PCTTU| have 
analyzed data from measurements of the decay rates of ^^Si and '^^Cl acquired 
at the Brookhaven National Laboratory (BNL; [Xlburger, Harbottle, and Norton| 
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(|1986p ) and from measurements of the decay rate of ^^^Ra acquired at the 
Physikahsch-Techniche Bundesanstah (PTB; |Siegert, Schrader, and Schoetzig|(|199g)) ).| 
finding evidence of smah but statisticahy significant annual variations in those 
rates. These findings led those authors to suggest that nuclear decay rates may be 
influenced (via some unknown mechanism) by the Sun. These analyses have been 
criticized by [Cooper] (gnOSl) , INorman et aL\ (POT]) , and lSemkow et al.\ (^005)) , but 
detailed responses to these criticisms have been presented by |Jenkins, Mundy, and Fischbach^ 
(pgrn]) and bv lJavorsekll et g/.K^UTU)) . 

The possibility of a solar influence on decay-rate variations has led us to search 
for periodicities that might be indicative of solar rotation. Our recent analyses 
(jSturrock et al.[ I2010al I2010bp of BNL and PTB data yield evidence in both 
datasets of a periodicity with a frequency in the neighborhood of 11.25 year~^[3 
The standard shuffle test ( [Bahcall and Press, 1991[ ) indicates that there is a 
probability of only lO^^*" of finding by chance both peaks at the same frequency 
in a search band of 10-15 year~^. We have also examined the data by means of 
a new test, called the "shake" test ( [Sturrock et ai, 2010b[ ), which yields results 
very similar to those obtained by the shuffle test. 

This periodicity is suggestive of the influence of a core or the environment 
of a core that has a sidereal rotation frequency of about 12.25 year~^. [The 
synodic rotation frequency — the frequency as observed from Earth — is lower 
than the sidereal (absolute) rotation frequency by 1 year~^.] The thin region 
at normalized radius approximately 0.7 that separates the radiative zone from 
the convection zone is known as the "tachocline," since it is the location of a 
sharp change in rotation rate. If there is a similar sharp change in rotation rate 
between the core, with sidereal rotation frequency at or close to 12.25 year~^, 
and the radiative zone, with a mean sidereal rotation rate of about 13.7 year~^ 
( [Schou and et al., 1998[ ), there would be a second or "inner" tachocline, which 
could conceivably have properties similar to those of the "outer" tachocline. 

There is supporting evidence for this conjecture. A combined analysis (|Sturrock[ [ 
I200 8','2009l) of low-energy solar neutrino data from the Homestake (jCleveland et a/.[ | 
M9 8. Davisl UM^ and GALLEX (jAnselmann et aU [TOMl UM^ experiments 
and of total solar irradiance data from the ACRIM experiment (|Willson[ 119791 
I200ip yields evidence of a periodicity of approximately 11.85 year~^, corre- 
sponding to a sidereal rotation frequency of 12.85 year~^. Also, a recent anal- 
ysis of Mount Wilson solar diameter measurements yields evidence of oscilla- 
tions originating in a region with sidereal rotation frequency of 12.08 year~^ 
dSturrock and Bertello, 2010] ). These inferred rotation frequencies are not the 
same, but they are both lower than current estimates of the rotation rate of the 
radiative zone ( [Schou and et al., 1998[ ), so that they could both originate in a 
tachocline. 

The well-known solar Rieger periodicity (which was originally discovered in 
gamma-ray flare data; [Rieger et al.\ I1984p has a period of about 154 days, cor- 
responding to a frequency of approximately 2.37 year~^. We have proposed that 



^An earlier anal ysis ||Sturrock et a l, 2010b^ led to the estimate 11.93 year~^ for BNL, but our 
recent analysis jSturrock et al., 2010aj shows that this estimate arose from the convolution of 
two separate peaks in the power spectrum, one at 11.24 year~^ and the other at 13.08 year~^. 
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this periodicity may be attributed to an r-mode oscillation that occurs in the 
known tachocline ( [Sturrock et ai, 1999[ [Sturrock, Caldwell, and Scargle, 20"06{ 



Sturrock, 2009), an interpretation that receives strong support from the recent 



article by Zaqarashvili et a/.] (|2010p . 



These considerations suggest that we may test the hypothesis of a solar in- 
fluence on nuclear decay rates as follows. If there is indeed an inner tachocline, 
it may be the seat of a Rieger-like oscillation at a frequency determined by the 
band of rotation frequencies of that region. This leads to a specific proposal that 
we present in Section 2. Our search for this periodicity is presented in Section 
3. We discuss the results in Section 4. 



2. Proposal 



r-mode oscillations (Papaloizou and Pringle, 1978 Provost, Berthomieu, and Rocca, 1981H 
|Saio, 1982] Wolff and Blizard, 1986[ ) are retrograde waves in a rotating fluid with 



frequencies determined, to good approximation, in terms of I and m (two of 
the three spherical harmonic indices) and the absolute or "sidereal" rotation 
frequency, . The values of I and m are restricted by 



1 = 2,3,..., m = l,2,...,l (1) 

The r-mode frequencies, as they would be measured by an observer co-rotating 
with the Sun, are given to good approximation by 

Equation (2) yields the 2.37 year~^ Ricger frequency for the combination 
Z = 3, m = 1 and vn = 14.22yr^^, while other values of l,m and vj^ and are 
disfavored. We find, from estimates of the internal solar rotation rate obtained 



from helioseismology (Schou and et al., 1998), that the Sun has this rotation 



rate at normalized radius 0.71, placing it in the tachocline. 

Taking into account (a) the inferred sidereal rotation rate of 12.25 ycar~^ with 
an uncertainty of 0.13 year~^ derived from BNL and PTB data, (b) the estimated 
sidereal rotation rate of 12.85 year~^ derived from solar neutrino and total solar 
irradiance data, and (c) the estimated sidereal rate of 12.08 year"-'^ derived from 
Mount Wilson diameter data, we adopt 12.0 year^^ for the lower limit of the 
sidereal rotation frequency of the inner tachocline. For the upper limit, we adopt 
the mean rotation rate of the radiative zone, 13.7 year~^. Adopting the same 
values of 1 and m {I = 3, m = 1) as in our interpretation of the Rieger oscillation, 
we see from Equation (2) that the corresponding search band for the frequency 
of a Rieger-typc periodicity associated with the hypothetical inner tachocline is 
2.00-2.28 year-i. 

These considerations lead to the following proposal: The hypothesis of a solar 
influence on nuclear decay rates leads us to expect possible periodicities in BNL 
and PTB data in the range 2.00-2.28 year~^ . 
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Figure 1. Power spectrum formed from the BNL data. As mentioned in the text, this figure 
incorporates additional data obtained between 1986 and 1990, which were not utilized by 
Jenkins ct al. (2009) or by Fischbach et al. (2009). The arrows indicate the search band 2.02 
year-i to 2.28 year"!. The peak is found at 2.107 ± 0.004 year^i with power S = 10.09. 



3. Analysis 



Wc show in Figures 1 and 2 sections of the power speetra recently derived 
from BNL and PTB datasets, respectively (jSturrock et aZ.i I2010all2010bp . The 
complete BNL dataset comprises 366 measurements in the time interval 1982.11 
to 1989.93. We point out that, due to a change in the apparatus in 1986, 



Alburger, Harbottle, and Norton] (|T986| utilized only data from 1982 to 1986 
in their determination of the ^^Si half-life. However, we have chosen to analyze 
the complete dataset in our power spectrum analysis, since the slight disconti- 
nuity in the dataset should not affect the frequency content (although it would 
affect the determination of the "^^Si half-life, which was the objective of the BNL 
experiment). 

Following the advice of the BNL experimenters, we have examined the ratio 
of the measurements of the decay rates of 32Si and 36C1, each detrended to take 
account of the mean decay rate. We also examine the 1,966 PTB measurements 
of the decay rate of 226Ra, similarly detrended. Power spectra were formed 



using a likelihood procedure ( Sturrock et al, 2005a I that is equivalent to the 
Lomb-Scargle procedure ( Lomb, 1976 Scargle, 1982 1. The power spectra are 
complex, showing strong features at 1 year~^ and at other frequencies (such 
as the neighborhood of 1.7 year~^ in Figure 2). We plan to discuss these low- 
frequency features in a future article, but we here focus on features within the 
search band. 

The BNL power spectrum shows a peak at 2.11 year~^ with power SI = 10.09. 
In order to obtain an error estimate, we adopt the convention of determining the 
"1-sigma" values. For a normal distribution of a quantity, the 1-sigma error is the 
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Figure 2. Power spectrum formed from the PTB data. The arrows indicate the search band 
2.02 year~^ to 2.28 year~^. The peak is found at with power S = 25.83. 



displacement that leads to an increase in the probability of the null hypothesis by 
the factor e°'^. In power-spectrum analysis, the probability of the null hypothesis 
(that the peak is due to normally distributed random noise) is given by 
( [Scargle, 1982 1. Hence the 1-sigma offset is that at which S is decreased by 0.5. 
In order to determine this offset for the BNL power spectrum, it is necessary 
to use a frequency resolution of only 0.001 ycar^^. We then obtain the estimate 
2.108 ± 0.012 year^^ for the location of the peak. 

A similar analysis of the PTB data yields a peak at 2.11 year^^ with power 
S2 = 25.83. Adopting a frequency resolution of 0.001 year"-'^, we obtain the 
estimate 2.107 ± 0.005 year^^ for the location of the peak. Taking into account 
the error estimates, the frequencies of the BNL and PTB peaks are clearly 
compatible. 

We may examine the correlation between these two power spectra by forming 



the joint power statistic ( Sturrock et al., 2005b ) given, to good approximation, 

by 



J = 



1.943X2 



(0.650 + X) 
where X is the geometric mean oi Si and 6*2: 



(3) 



X = {8182^' (4) 

This statistic has the property that if 5*1 and 5*2 are each distributed ex- 
ponentially (as is the case for a time series derived from normally distributed 
random noise; [Scargle[ [TOS^ . then J also is distributed exponentially. We show 
J as a function of frequency in Figure 3. The peak at 2.11 year~^ has the 
value J = 30.65. On taking account of our recent discussion of the false-alarm 
probability of peaks in power spectra ( [Sturrock and Bertello, 2010] ) , we estimate 
that probability to be 1.3 x 10"^'^. On adopting the finer frequency resolution, 
we find that the frequency of the peak value of the JPS is 2.107 ± 0.005 year~^. 
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Figure 3. The joint power statistic formed from the BNL and PTB power spectra. The arrows 
indicate the search band 2.02 year~^ to 2.28 year^^. The peak is found at 2.107±0.005 year~^ 
with joint power statistic J = 30.65. 
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Figure 4. Histogram formed from 10,000 shuffle simulations of the joint power statistic formed 
from the BNL and PTB power spectra. The vertical red hne indicates the value (30.65) of the 
joint power statistic for the peak at 2.11 year~^. 



Although the separation between peaks is given approximately by l/T, where 
T is the duration of the dataset, we have seen that the precision of an individual 
peak depends critically on the peak power, and can be much less than the l/T 
value. The frequency estimates derived from BNL data, from PTB data, and from 
the joint power statistic all differ significantly from 2.00 year~^, the frequency 
of the first harmonic of the annual oscillation. 
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Figure 5. Logarithmic display of the results of 10,000 shuffle simulations of the joint power 
statistic formed from the BNL and PTB power spectra. The projection of this curve indicates 
that the probability of obtaining by chance a value of J as large as the actual value (30.65) is 
about 10-12- 



The significance of a peak in the power spectrum and the properties of the 
joint power statistic are based on an assumed exponential form for the distribu- 
tion function of the power, which is appropriate if the non-signal measurements 
(the "noise") conform to a normal distribution (Scargle, 1982). We examine the 
validity of this assumption for the BNL and PTB datasets in the Appendix, 
where we show that mapping the measurements onto a strictly normal distribu- 
tion has a negligible effect on the power spectrum, from which wc may infer that 
any departure of the measurements from a normal distribution has no significant 
adverse effect on our analysis. 

However, in order to obtain a robust significance estimate, we have applied 
the shuffle test ( [Bahcall and Press, 1991[ ) that does not depend on any assumed 
form of the power distribution function. We have computed J for 10,000 Monte 
Carlo simulations generated by the shuffle procedure, shuffling both datasets 
and noting the maximum value of J in the frequency band 2.00-2.28 year~^ for 
each simulation. The results from the shuffle test are shown in histogram form 
in Figure 4 and in a logarithmic display in Figure 5. A projection of the results 
shown in Figure 5 indicates that the probability of obtaining by chance a value of 
J as large as the actual value (30.65) is about 10~^^, slightly more conservative 
than the false-alarm probability. 

We have also carried out simulations using the shake test ( [Sturrock et ai, 2010b[ ),| 
which has the advantage over the shuffle test that it is less sensitive to low- 
frequency modulations and secular trends. However, the results of the shake and 
shuffle tests prove to be virtually identical, indicating that there is negligible 
probability of obtaining by chance a value of the JPS as large as or larger than 
the actual value. Wc find, therefore, that this analysis is consistent with our 
proposal. 
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As a further precaution against the possibihty that the osciUations at 2.11 
year~^ represent the first harmonic of the strong annual osciUations, we have 
carried out the foUowing test. We have formed a "template" of the actual an- 
nual variation of the BNL data, and divided the actual measurements by the 
corresponding values of the template. This procedure should remove all effects 
of the annual oscillation, including harmonics. We find that the peak at 2.11 
year~^ remains in the power spectrum. We have applied the same procedure to 
the PTB dataset with the same result. This test indicates that the peaks at 2.11 
year"^ in the BNL and PTB power spectra may not be attributed to a harmonic 
of the annual oscillation. 



4. Discussion 

It would appear that this study further strengthens the case for a solar infiu- 
ence on nuclear decay rates. It does not offer further insight into the physical 
mechanism of this influence, but it supports the inference that the origin of the 
influence is to be found below the radiative zone — either in or near the core. 

The results of our analysis are therefore relevant to the question of the rotation 
rate of the solar core. It has to date been difficult to determine this quantity by 
means of helioseismology, and the results have been conflicting: |Chaplin et al.\ 
P999p reported evidence for a downturn in the rotation rate below normalized 
radius 0.15, but [Chaplin et ar](|200ip later estimated the rate to be close to that 
of the radiative zone. On the other hand, IGarca et al.\ (|2007p have estimated 
that the rate may be three to five times faster than that of the radiative zone. 
IHowel (|2009p , in her recent review of solar rotation, concluded that "the rotation 
of the inner core remain(s) unclear." The results presented in this article provide 
further evidence (in addition to ISturrock et ail ISturrock et al.\ (|2010al I2010bp ) 
that rotation of the solar core is slower than that of the radiative zone. However, 
we are in a very early stage of learning about the structure and rotation of the 
deep solar interior. Any model that we can propose at this time will no doubt 
need to be revised in the future. 

One may ask why r-mode oscillations should be preferentially excited in a 
tachocline. |Zaqarashvili et aL|(pUTU|) . who also attribute the Rieger periodicity 
to an r-modc oscillation in the outer tachocline, have shown that r-modes become 
unstable in a layer where there is both a toroidal magnetic field and latitudinal 
differential rotation, and that this instability favors m — 1 modes. However, 
this interpretation does not explain why the instability should be confined to 
the tachocline, since similar conditions may be presumed to exist throughout 
the convection zone. Our current conjecture is that r-modes are unstable in the 
presence of radial differential rotation. We plan to investigate this hypothesis 
in the near future. Wc plan also to search for evidence of other r-modes, in 
addition to the mode with / = 3,to = 1, which offers an explanation of the 
Rieger oscillation. 
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Appendix 



A. The Rank-Order Normalization (RONO) Procedure 

The power-spectrum analyses in this article have been carried out using a likeli- 
hood procedure (Sturrock et ai, 2005a[) that is equivalent to the Lomb-Scargle 



procedure ( Lomb, 1976 Scargle, 1982[ ). However, these procedures arc based on 



the assumption that the background (non-signal) data have a normal distribu- 
tion. Since there is no reason to believe that the background data for the BNL 
and PTB experiments are distributed strictly in this way, we here examine the 
question of whether the departures of the BNL and PTB datasets from normal 
distributions are of any consequence. We address this question by the following 
procedure, which we refer to as the "RONO," for "Rank-Order Normalization," 
procedure. 

The dataset comprises a sequence of measurements, Xn,n = I, . . . ^N, which 
we denote by {x}, taken at times tn,n ~ 1,. . . ,N, which we denote by {t}. We 
now arrange the sequence of measurements in ascending order. We denote the 
re-arranged measurement sequence by {x} = R{x}. 



f = erf ig),g = erfinv (/) (5) 

We now define {y} by 

y = erfinv (/„) (6) 

where 



fn = ,n = l,...,Af. (7) 

iV + 1 

We have applied this transformation to the BNL data, and then applied the 
likelihood power spectrum analysis to the normalized data {y}. The result is 
shown in Figure 6, in which the original power spectrum (derived from the 
measurements {x}) is shown in blue, and the new power spectrum (derived from 
the normalized measurements {y}) is shown in red. The normalized data lead to 
a peak at 2.11 year~^ with peak power 10.43. We see that the frequency of the 
peak is indistinguishable from the value in the original analysis, and the peak 
power is changed only slightly (from 10.09). 

The same operation on the PTB data leads to a power spectrum that has a 
peak at 2.11 year"-'^ with power 25.50 — again no change in frequency, and only 
a slight change in power from the original estimate of 25.83. 
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Figure 6. Tiic power spectrum formed from BNL data, after the application of the RONO 
procedure, is shown in red. The original power spectrum is shown, for comparison purposes, 
in blue. The original peak had power 10.09 at frequency 2.11 year~^. The power spectrum 
derived from the normalized data has power 10.43 at the same frequency 2.11 year"^. 
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